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Abstract 

We review recent progress in implementing high-precision chiral two- and three-body forces in 
nuclear many-body systems beyond light nuclei. We begin with applications to finite nuclei, which 



we study through the nuclear shell model and self-consistent mean field theory. We then turn our 
attention to infinite nuclear matter treated within the framework of Landau's theory of normal 
Fermi liquids. 

2 1 Introduction 

S3 

In recent years, the application of chiral effective field theory methods to low-energy nuclear dynamics 
has resulted in the development of high-precision two- and many-body forces consistent with the sym- 
metry constraints imposed by the fundamental theory of strong interactions, quantum chromodynamics. 
Q\ ', For recent reviews on this subject, see refs. [HE]- Although such interactions have been highly success- 
^ \ ful in ab initio calculations of light-nuclei scattering observables, binding energies and spectra [3], HIE], 
there remains much work to be done to implement three-nucleon forces in studies of medium-mass and 
heavy nuclei. In this contribution we discuss our recent work to understand the properties of finite 
nuclei and infinite nuclear matter from many-body perturbation theory employing realistic N 3 LO chiral 
two-nucleon forces supplemented with the leading-order N 2 LO three-body force. For this purpose it is 
.£h ! convenient to employ as well low-momentum nucleon-nucleon (NN) interactions [21 E], which exhibit 
^ ' favorable convergence properties in perturbative calculations [S]. 



2 In-medium effective nucleon-nucleon interactions 



To facilitate the implementaion of three-body forces in nuclear many-body calculations, it is useful to 
construct medium- dependent NN interactions that reflect the most important physical features of the 
underlying three-nucleon force. Here we focus on a background medium of isospin-symmetric nuclear 
matter and compute the density-dependent NN interaction to one-loop order from the N 2 LO chiral 
three-nucleon interaction. In fact, the effective interaction in nuclear matter with isospin asymmetries 
up to 5 np = (p n — Pp)/p — 0.2 is already well-approximated by the effective interaction for S np = 
[9]. For detailed discussions regarding the in- medium interaction in pure neutron matter, we refer the 
reader to refs. P \10\ . 

The leading-order chiral three-nucleon force includes three separate components of varying range. 
The two-pion exchange contribution V^ffl consists of terms proportional to the low-energy constants 
ci )3j 4, which arise already in the two-nucleon interaction and can therefore be constrained by fits to 
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Figure 1: In-medium NN interaction generated by the three different components of the chiral three- 
nucleon force. The short double-line symbolizes the filled Fermi sea of nucleons, and reflected diagrams 
are not shown. 

NN scattering phase shifts. The one-pion exchange component is proportional to the low-energy 

constant cd, while the short-range contact contribution vffl introduces the low-energy constant eg. 
Both cd and ce are commonly fit to the properties of A = 3, 4 systems. In the present work we use two 
different parameterizations of the low-energy constants ce, cd and 01,3,4: 

c B = -0.205, c D = -0.20, d = -0.81 GeV" 1 , c 3 = -3.2 GeV" 1 , c 4 = 5.4GeV~\ (1) 

c E = -0.625, c D = -2.06, c x = -0.76 GeV" 1 , c 3 = -4.78 GeV^ 1 , c 4 = 3.96 GeV" 1 , (2) 

which are consistent with the Idaho N 3 LO chiral NN interaction [TTJ and the low-momentum NN 
interaction V[ ow _k evolved to a decimation scale of A = 2.1 fm _1 |8], respectively. 

At one-loop order there are six distinct diagrams, shown in Fig. [TJ contributing to the in-medium 
NN interaction Vjy^r. The double-line symbolizes summation over the filled Fermi sea of nucleons, or 
equivalently, the medium insertion —2ix5(ko) 9(kf — \k |) in the in-medium nucleon propagator. For a 
background medium at rest, the in- medium scattering amplitude for two nucleons in their center of 
mass frame, N\(p) + N2(—p) — > N\(p + q) + A^— p — q), has the same general form as that given in 
free-space NN scattering. For a background medium boosted with respect to the center of mass frame, 
the effective interaction will in general depend also on the boost momentum. 

In Fig. [2] we plot the 1 5'o and 3 Si diagonal partial- wave matrix elements of the low-momentum NN 
potential V[ ow _k and the modifications resulting from the six components of evaluated at nuclear 
matter saturation density (po — 0.16 fm -3 ) with low-energy constants given in eq. (TS]). In general, the 
three contributions from can be a significant fraction of the free-space matrix elements of V[ ow _k 
and are considerably stronger than those arising from the one-pion exchange and contact three-nucleon 
forces. The pion self-energy correction (diagram (1) in Fig. [1]) gives rise to an enhancement of the bare 
l7r-exchange, which can be interpreted in terms of a reduced pion decay constant, f* 2 s — f% + 2c 3 p. 

In contrast, the vertex correction (diagram (2) in Fig. [T]) from V^j reduces one-pion exchange, an 
effect that can be approximated by introducing a reduced nucleon axial-vector constant g* A . Taken 
together, these two terms largely cancel in all partial waves, which leaves the Pauli-blocked two-pion 
exchange process (diagram (3) in Fig. [T| as a dominant contribution. The contributions from the 
one-pion exchange three-body force are quite small. Since cd is negative, the vertex correction 

(diagram (4) in Fig. [TJ) reduces the bare l7r-exchange, but only by about 16% at normal nuclear matter 
density. The pion-loop correction to the NN contact interaction (diagram (5) in Fig. [T]) acts only in 
the two relative S- waves as well as the S — D mixing matrix element. In relative 5*- waves, it reduces 
the repulsive contribution from the pure contact term (diagram (6) in Fig. [Q by about one-half. In 
the right plot of Fig. [2] we sum all six contributions, which together give rise to significant additional 
repulsion that increases with the density and which provides a mechanism for nuclear matter saturation 
when employing low- momentum potentials [S]. 

The resulting in-medium effective interaction has been employed already in shell model studies of 
finite nuclei, where it was found to play a key role in enhancing the /3-decay lifetime of 14 C [12] toward 
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Figure 2: Left figure: contributions to the S'-wave effective interaction in symmetric nuclear matter 
(at saturation density, po — 0.16 fm~ 3 ) resulting from the six diagrams shown in Fig. [TJ The solid 
black line shows the matrix elements of the free-space NN potential V[ ow _k, and the labels 1-6 denote 
the contributions from the corresponding diagrams in Fig. [IJ Right figure: contributions from all 
components of the in-medium NN interaction to the L = partial waves at p = p ,p Q /2. 



its known achaeologically-long value of 5730 years. An alternative approach for finite nuclei calculations 
is to evaluate the three-nucleon force in a shell model basis and subsequently sum over the filled orbitals 
in the closed nuclear core. Such a method has been used to explore the effects of three-nucleon forces 
in neutron-rich oxygen and calcium isotopes [13]. 

3 Microscopic energy density functionals 

To study the binding energies and charge radii of nuclei across the periodic chart, the most successful 
theoretical framework has been self-consistent mean field theory based on phenomenological nuclear 
energy density functionals [14|. A microscopic foundation for such functionals is provided within the 
framework of many-body perturbation theory by the density matrix expansion of Negele and Vautherin 
[15], in which the highly nonlocal expression for the energy is expanded in terms of local densities 
and currents (and their gradients) resulting in a generalized Skyrme functional with density- dependent 
couplings. 

Recently Gebremariam, Duguet and Bogner [16J have developed an improved density matrix ex- 
pansion that gives rise to a much improved description of the spin density matrix. The form of the 
density matrix they find (neglecting isovector terms) can be expanded in relative and center-of-mass 
coordinates, a and r, as follows: 



J2*«(r-3/2)*Ur + a/2) = *g- j x {ak f ) - ±- j x {ak f ) 

a aK f AK ] 



r - \pk) - \v 2 p 



3i 

+ ^rh(ak f )a- (ax J) + ... , (3) 

where the spherical Bessel function ji(x) = (sinx — xcosx)/x 2 . The quantities appearing on the 
right hand side of eq. ([3]) are the (local) nucleon density p(r) = 2kj-(r) /3n 2 = J2 a ^ai^^ai^), the 
(local) kinetic energy density r(r) = EaVf^f) • V\l/ Q (r) and the (local) spin-orbit density J{r) = 
«E a t(^ x Vf tt (r). 
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Figure 3: Contributions to the energy per particle E(p) of isospin-symmetric nuclear matter and to the 
spin-orbit strength parameter F so as a function of the nuclear density p. 



The energy density functional resulting from this improved density matrix expansion has been 
calculated to three-loop order including l7r-exchange, iterated l7r-exchange and irreducible 27r-exchange 
with intermediate A-isobars in ref. [17] and with a microscopic N 2 LO nucleon-nucleon potential in ref. 
[T8] . In our work [19], we have employed a microscopic N 3 LO chiral NN potential supplemented with 
the N 2 LO chiral three-nucleon force. The two-body interaction is comprised of long-range one- and two- 
pion exchange contributions together with a set of contact terms contributing up to fourth power in 
momenta. The potential, referred to as the N 3 LOW chiral NN interaction [20], has a sharp momentum- 
space cutoff A = 2.1fm _1 corresponding to the scale at which low-momentum interactions become 
universal and strongly perturbative. The low-energy constants eg, cp and ci j3j 4 of the leading-order 
chiral three-nucleon interaction are those of eq. 

Up to second order in spatial gradients, the energy density functional relevant for N = Z even-even 
nuclei reads: 



£\p,r, J]=pE(p) + 



3 _ n -ir 1 k 2 



f 



2M 4M 3 



+ F T (p) 



+ (Vp) 2 F v (p) + Vp ■ JF so {p) + J 2 Fj(p) . (4) 



Due to the nonlocalities present in the chiral contact interactions, it is convenient to compute the 
density- dependent strength functions E(p), F T (p), F d (p), F so (p) and Fj(p) separately for the finite- 
range pion-exchange and the zero-range contact components in the chiral N 3 LO NN interaction. The 
leading-order three-nucleon force is computed under the assumption that the relevant product of density- 
matrices can be represented in momentum space in a factorized form: r(p*i, <fi) T(p 2 , q 2 ) r(p~3, —q\ — q\). 

In Figs. |3H5]we plot the density dependence of the five strength functions E(p), F T (p), Fd(p), F so (p) 
and Fj(p). At first order in perturbation theory, low-momentum interactions are typically underbound 
after including the leading-order chiral three-nucleon force [8], and this is reflected in the left plot of 
Fig. [3J where it is seen that both the saturation density and binding energy are too small. In the right 
plot of Fig. [3] we show the spin-orbit strength, which receives a very large positive contribution from the 
N 2 LO chiral three-body force that increases the strength at p — po/2 beyond values typical of Skyrme 
functionals. The strength function F T (p) (shown in the left plot of Fig. H|) is related to the effective 
nucleon mass M*, shown in the right plot of Fig. [H At saturation density, the effective nucleon mass 
is M* ~ 0.7 Mjy, within the range of values used in phenomeno logical functionals. Finally, we consider 
Fv(p) and Fj(p) shown in Fig. The former encodes the energy due to density gradients at the Fermi 
surface, and at half nuclear matter saturation density it achieves values slightly below those extracted 
from Skyrme functionals. The strong variation at low densities of Fj(p) originates from the dominant 
l7r-exchange contribution [T7J. The resulting microscopic energy density functional is an encouraging 
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Figure 4: Contributions to the strength parameter F T and effective mass M* as a function of the nuclear 
density p. 




Figure 5: Contributions to the strength parameters Fy and Fj as a function of the nuclear density p. 



start and leaves room for second-order perturbative contributions, which are expected to improve in 
particular both the energy per particle E(p) and the spin-orbit strength F so [T9] . 



4 Quasiparticle interaction in nuclear matter 

Landau's theory of normal Fermi liquids [2T] can be employed to study the bulk equilibrium and trans- 
port properties of strongly-interacting normal Fermi systems at low temperatures. The theory describes 
the low-energy excitations above the interacting ground state in terms of quasiparticles, which retain 
certain features of noninteracting particles and in a sense are weakly interacting. The quasiparti- 
cle interaction encodes bulk properties of the medium and dynamical properties of the quasiparticles 
themselves. It can be computed in many-body perturbation theory by functionally differentiating the 
ground-state energy twice with respect to the quasiparticle distribution function and has the form 

1 oo 

HPiiPi) = irr £ [ p L + ■ r 2 + (G L + G' L n ■ 7^1 ■ 3 2] P L (cas0), (5) 

7V L=0 

where we have expanded the angular dependence in Legendre polynomials, factored out the density of 
states iVo at the Fermi surface and included only the central components of the quasiparticle interaction. 
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Table 1: Dimensionless Fermi liquid parameters at nuclear matter saturation density (corresponding to 
hp = 1.33 fm _1 ) obtained from chiral two- and three-nucleon interactions as described in the text. 



Recently, we have employed chiral two- and three-body forces in a systematic study of the quasi- 
particle interaction in isospin-symmetric nuclear matter [22| 123] . Already a number of observables (see 
eq. OH]) for the relation to Fermi liquid parameters) were well described in the second-order perturbative 
calculation that included two-nucleon forces only [22]. In particular, the values of the quasiparticle 
effective mass M*, the isosp in- asymmetry energy (3 and the spin-isospin response parameter g' NN were 
found to be M* /Mjy — 1.0 — 1.1, ft 31 — 32 MeV and g' NN ~ 0.65 — 0.75, all within empirical ranges. 
However, both the compression modulus /C ~ —200 MeV and the anomalous orbital ^-factor Sgi ~ 0.1 
were found to differ appreciably from the empirical values /C ~ 200 — 300 MeV and 5gi ~ 0.20 — 0.26 
inferred from studies of giant resonances in heavy nuclei. 

As a first step toward the consistent inclusion of the N 2 LO chiral three-nucleon force, we have com- 
puted the first-order perturbative contribution to the L — 0, 1 Landau parameters of the quasiparticle 
interaction. Two sets of low-energy constants, given in eqs. ([I])-©, were employed, which allowed for 
an estimate of the theoretical uncertainty at this order in the calculation. The diagrammatic contri- 
butions to the quasiparticle interaction are shown in Fig. HJ but since the two quasiparticles lie on the 
Fermi surface \pi\ = \p2\ = kp, the kinematics in the present case are different than those considered 
in Section [2j We show in Table [1] the results from the second-order calculation including two-nucleon 
forces only (labeled "2iV" ) and the results obtained by including as well the one-loop contributions from 
the N 2 LO chiral three-body force (labeled "3iV"). We notice that the most visible effect is a dramatic 
increase in the isotropic spin- and isosp in- independent parameter F , which was large and negative 
(leading to a negative compression modulus K) when computed with only chiral and low-momentum 
NN interactions, but which now attains values in reasonable agreement with empirical constraints from 
giant monopole resonances. In comparison, the contributions to the other Landau parameters from 
the chiral three-nucleon force considered here are much smaller, and most observables remain in good 
agreement with their empirical values: 

M* 

Effective mass : — - = 1 + FJ3 = 0.98 ± 0.05, [0.7 - 1.0] 

M N 

Anomalous orbital ^-factor: Sg t = — ^ ^- = 0.07 ± 0.02, [0.20 - 0.26] 

y y 6(l + Fi/3) ' L J 

Compression modulus : K = — J (1 + F ) = (370 ± 160) MeV, [200 - 300] MeV 

k 2 

Isospin asymmetry energy : /3 = ^-(1 + F^) = (30 ±1) MeV, [30 - 36] MeV 

AM 2 

Spin-isospin response : g' NN = — — —G' = 0.55 ± 0.03 [0.6 — 0.7]. (6) 



o 



In these equations, g n N ~ 13.2 is the strong irN coupling constant and the values in brackets represent 
empirical estimates (for further details, see ref. [23]). Moreover, the theoretical uncertainty resulting 
from different choices of cutoff scale and low-energy constants is reduced with the inclusion of three- 
nucleon forces. The one observable which remains in disagreement with its empirical value is the 
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anomalous orbital g-factor 5gi, which lies well below the value 5gi = 0.20 — 0.26 extracted from a sum- 
rule analysis of giant dipole resonances. As suggested in ref. [22], this could be remedied if higher-order 
perturbative calculations reduce even moderately the quasiparticle effective mass M*. Finally, we notice 
that all L = 1 Landau parameters decrease with the addition of three-nucleon forces, which results in 
an effective interaction of apparent short range. 



5 Conclusion 

In this talk we have focused on our recent efforts to incorporate the leading-order chiral three-nucleon 
force in many-body perturbation theory calculations of dense nuclear systems. This work has set the 
foundation for future shell model and mean field theory calculations of finite nuclei, and our investiga- 
tions of infinite isospin-symmetric nuclear matter can be broadened to include pure neutron matter as 
a first approximation to neutron-star matter. In the future we look forward to employing the recently 
completed subleading chiral three-nucleon force [23], which will allow for consistent calculations up to 
order N 3 LO in the chiral expansion. 
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